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Abstract 

We develop a complete Dirac-Hartree-Fock-Bogoliubov approximation to the 
ground state wave function and energy of finite nuclei. We apply it to spin- 
zero proton-proton and neutron-neutron pairing within the Dirac-Hartree- 
Bogoliubov approximation (we neglect the Fock term), using a zero-range 
approximation to the relativistic pairing tensor. We study the effects of the 
pairing on the properties of the even-even nuclei of the isotopic chains of Ca, 
Ni and Sn (spherical) and Kr and Sr (deformed), as well as the N=28 isotonic 
chain, and compare our results with experimental data and with other recent 
calculations. 
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I. INTRODUCTION 



The advances made at radioactive nuclear beam facilities provide us with an increasing 
amount of information on nuclei far away from the stability line. We now have access to 
experimental measurements of the masses, radii and deformations of unstable nuclei in a 
wider region of the nuclear chart than ever before. Studies of exotic nuclei have revealed 
new features, such as neutron halos [0 and neutron skins [§,[3| and bring new perspectives 
to nuclear physics 0. Planned facilities around the world plan to study unstable nuclei up 
to the r-process region or even beyond. This will enable us to investigate where and how 
exotic phenomena of nuclear structure appear in the region far from the stability line . 

In recent years relativistic many-body theories have been applied to nuclei and nuclear 
matter with remarkable success [0^01 • The relativistic Brueckner-Hartree-Fock (RBHF) 
theory has been shown capable of reproducing the saturation properties of nuclear matter 
using interaction parameters that describe the two-nucleon bound state and scattering |TU 



Its phenomenological version, the relativistic mean field model (RMF), has been successfully 
applied to the description of many of the ground-state properties of stable and unstable nuclei 
|TT^p4| and has been shown capable of simultaneously reproducing the properties of stable 
and unstable nuclei over a wide mass range of the periodic table . 



The relativistic mean field approach has also been applied to describe the structure 
of very exotic nuclei. In these calculations, the pairing interaction has been neglected or 
simply treated by a nonrelativistic BCS type of approximation |2^J2^]. For nuclei on or 
near the stability line, the BCS approach provides a reasonably good description of pairing 
properties. However, for drip-line nuclei, the Fermi level is close to the continuum and the 
coupling between bound and continuum states should be taken into account explicitly. The 
pairing correlations and the mean fields should be calculated simultaneously in order to 
obtain an accurate description of the ground state properties of drip-line nuclei. 

A more precise relativistic description of pairing correlations, a Dirac-Hartree-Fock- 
Bogoliubov (DHFB) approximation, was developed some time ago |27] and, more recently, 
applied to nuclear matter calculations [E28H3T|. The first calculations in nuclear matter fur- 
nished pairing gaps much larger than those obtained in nonrelativistic calculations using 
realistic interactions. |28| The results suggested that the meson-exchange interactions ad- 
justed to describe nuclear matter saturation might not be adequate for the particle-particle 
pairing channel. However, these calculations did not perform a complete self-consistent 
DHFB calculation. Rather, they combined a relativistic Hartree mean field calculation with 
a nonrelativistic calculation of the pairing field using the nonrelativistic reduction of the 
meson-exchange potential. Self-consistent calculations including the full Dirac structure of 
the self-energy and pairing fields resulted in more reasonable values for the pairing gap. 
29| , |30|| Yet these DHFB calculations, using various sets of meson-exchange parameters that 



furnish similar good descriptions of nuclear saturation, still resulted in values of the pairing 
gap that were not consistent with nonrelativistic calculations nor consistent among them- 
selves. In Ref. pTJ, these discrepancies were resolved by associating the pairing gap of 



each set of meson-exchange parameters with the set's description of low-energy two-nucleon 
scattering in the 1 So channel. When each of the parameter sets was supplemented with a 
momentum cutoff, so that they each described the two-nucleon 1 So virtual state correctly, 
they all provided consistent values for the pairing gap that very close to the nonrelativistic 
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ones, at low densities. Thus, the conclusion of Ref. |28| was found to be correct afterall. 
The meson-exchange interactions adjusted to describe nuclear matter saturation are indeed 
inadequate for the pairing channel, when they do not describe the two-nucleon 1 So channel 
well. In the case of nuclear matter, however, we have found that this inadequacy can be 
easily remedied. 

Even after the different sets of meson-exchange parameters have been augmented so 
as to describe the two-nucleon 1 So channel, they still differ slightly among themselves at 
densities greater than about one-eight of the saturation density, where mean field effects 
begin to play a role. Most importantly, they all differ in a consistent manner from the gap 
functions obtained using realistic nonrelativistic two-nucleon interactions. The relativistic 
pairing gap consistently falls to zero at lower density (by about one-half) of that for which it 
disappears in nonrelativistic calculations. Much as the Dirac structure of the nuclear mean 
field weakens the effective interaction as the density increases, leading to saturation, so the 
Dirac structure of the pairing field weakens the effective pairing interaction with density, 
leading to the suppression of the pairing field relative to the nonrelativistic one. An estimate 
of the relativistic effects on quasi-deuteron pairing shows that the same effects also suppress 
quasi-deuteron pairing at the saturation density |p2| . We thus expect a density-dependent 



suppression of the pairing field to be a general feature when we take the Dirac structure of 
the field into account. 

The relativistic Hartree mean field + nonrelativistic pairing field approximation of 



Ref. p8| can and has been extended to finite nuclei. This hybrid relativistic Hartree- 
Bogoliubov (RHB) approximation, using a non-relativistic finite-range Gogny interaction 
in the pairing channel, has been applied extensively, and quite successfully, to the calcula- 
tion of the ground state properties of spherical [|33H4T|1 and deformed nuclei [Q . We would 
expect reasonable results from such a model as long as the mean field is well described by 
the meson exchange interaction and the 1 So two-nucleon channel is well described by the 
effective pairing interaction. This is indeed the case in the references cited here. However, 
we might expect its pairing interaction to be less reliable than a meson-exchange one when 
extrapolated far beyond the stability valley, the condition for which its parameters have 
been adjusted. 

Aside from simply being more pleasing aesthetically, we are motivated to retain the 
Dirac structure of the pairing interaction in our calculations by the success of our studies of 
1 So pairing in nuclear matter and by their differences from the nonrelativistic results. The 
success of our nuclear matter calculations lead us to believe that a good description of nuclei 
can be obtained using a Dirac pairing interaction based on a meson-exchange one. Such a 
description would be of general interest, since the Dirac structure of a meson-exchange 
interaction could be extended beyond the valley of stability with better reliability than 
could a nonrelativistic interaction. A nonrelativistic effective interaction must incorporate 
the density dependence furnished by the Dirac structure of a meson-exchange interaction 
into its effective parameters. This density dependence changes as we consider nuclei further 
and further from the stability valley, thus limiting the range of application of a effective 
interaction adjusted to describe the conditions in the valley. 

The differences between our results and nonrelativistic ones in nuclear matter suggest 
that interesting effects of the Dirac structure of the pairing interaction might appear even in 
nuclei within the stability valley. The suppression of pairing in the saturated interior region, 
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in DHFB calculations, similar to that observed in nuclear matter, would result in greater 
localization of the pairing field on the nuclear surface than is the case in nonrelativistic 
or RHB calculations. This, in turn, would tend to make the nucleus more rigid and tend 
to diminish its deformation. We then would expect to observe a smaller pairing field in 
DHB calculations that furnish the same deformation as nonrelativistic or RHB ones, or a 
smaller deformation in DHB calculations that furnish the same value of the pairing field as 
nonrelativistic or RHB ones. 

Here we present a fully relativistic Dirac-Hartree-Bogoliubov (DHB) approximation for 
axially deformed nuclei. We use a direct extension of the DHFB approximation developed 
in Ref . |29| . We neglect the exchange term in the self-energy but retain the Dirac structure 



of the pairing interaction and pairing field. The neglect of the exchange term seems a 
reasonable approximation here, as we restrict our calculations to the exchange of a, p and 
lo mesons, for which the exchange effects can be fairly well accomodated by adjustments in 
the coupling constants. The characteristics of the pairing field and of its Dirac structure 
will form the center of our study. Our goal is to extend our successful DHFB description of 
nuclear matter to a DHB description of finite nuclei. 

We use a local (zero-range) approximation to the meson-exchange interaction in the 
pairing channel. A zero-range approximation to the pairing interaction can be justified 
by an analysis of the length scales in a typical DHB calculation. Taking oscillator wave 
functions as a guide, we can estimate the wavelength in one dimension of a wave function of 
quantum number n (in that dimension) as A n ~ 4Abo/\^n, where 60 is the oscillator length, 
bo ~ A 1 / 6 fm, and A is the mass number of the nucleus. In the fairly extreme case of the 
N=16 shell of 16 O, we find a wavelength of Ai6 ~ 1.75 fm. The wavelength is larger in 
heavier nuclei or in lower shells. On the other hand, the range of a typical meson-exchange 
interaction is r x = hc/m, where m is the meson mass. For the lightest of the mesons that 
we include, the a meson, we find a range of r a < 0.5 fm, substantially smaller than the 
wavelength Ai6- Even the term of longest range of the Gogny interaction has r 2 = 1.2 fm 
only. Thus we find that the interaction range is typically smaller than the length scale of 
the most energetic level, making the effects due to the finite range of the interaction small. 
This analysis has, in fact, been confirmed numerically by Meng | 38| , who found zero-range 
and finite-range Gogny RHB calculations to yield almost identical results. 

In Section I of the following, we derive the self-consistency equations for the self-energy 
and pairing fields. In Section II, we discuss the properties of static solutions to the equations 
and, in Section III, reduce these to the form in which we use them. In Section IV, we briefly 
describe the numerical method used in the calculations and, in Section V, we discuss the 
results of the calculations. We summarize and conclude in Section VI. 



II. THE MEAN FIELD EQUATIONS FOR S AND A 

In this section we present the Lagrangian of the model and the covariant equations for 
the self-energy and pairing mean fields, £ and A. 

We designate by ip(x), a(x), u^(x), p^^x), and A M (x), the field operators at the point x 
associated to the nucleons and mesons a, 00, p, and 7 respectively. The quantum numbers 
( ■/*", T) for each meson with spin J, intrinsic parity tt, and isospin T are 
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<x(0 + ,0), w(l-O), p(l",l), 7(1",-)- 

We designate the effective meson-nucleon coupling constants by g s , ^ and g p and the re- 
spective bare masses by m s , m v and m p . The interaction of the massless 7 does not conserve 
isospin, coupling to the protons alone with coupling constant e. (We will ignore the inter- 
action with the anomalous magnetic moments of the nucleons). The nucleon bare mass is 
M and we assume in the present model that the nucleons and mesons are point-like. We 
will not take the ir meson into account, although its effects could be important in the ex- 
change terms we will consider. We will assume that these effects can be described through 
the effective couplings of the mesons we do include. These assumptions are typical of the 
simplest meson-exchange models of nuclear structure. 
The Lagrangian density is given by 

£ — £0 + £int 1 
where C is the free Lagrangian density 

C (x) = ^(x)[ip - M]4>(x) + U p a(x)d»a{x) - U(a(x)) - -F pv F^ 

-In^Or + l -m 2 v ^{x)^{x) - \G, V ■ + \m%{x) ■ f?(x) , (1) 
with vector field tensors 

F 

1 (IV 
—* 

G pu 

and a nonlinear a potential, 

U(a(x)) = \mla{xf + l -g,a{xf + \g,a{xf . (2) 

The baryon spinor ip{x) has four Dirac components for each of two isospin projections - 
m t = 1/2 for protons and m t = —1/2 for neutrons - for a total of eight components. We 
have included the cubic and quartic terms of the scalar field a(x) in the free Lagrangian 
density as we will only consider their contributions to the scalar mean field. For this purpose, 
we may formally include them in the 'free' scalar meson propagator, when it is convenient 
to do so. 

We take the interaction terms in the Lagrangian density to have the simplest possible 
form consistent with their Lorentz and isospin structure, 

£int(x) = g s ip{x)a{x)ilj{x) - g v tl)(x)-i^uj' x (x)'il)(x) - ^ p ^(x)7 M f • p»(x)ijj(x) 

-e^-^pl^A^ix). (3) 
In particular, we will not consider tensor couplings of the vector mesons. 



— d p A v d u A fl 
= dpWv - d u u p , 
= du.Pi> ~ dupu , 
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We wish to characterize the average effect of the interactions of a nucleon with the other 
nucleons through an effective single-particle Lagrangian, L e // , given in terms of the two 
fields, S and A. The self-energy E describes the average interaction of a nucleon with the 
surrounding matter. The pairing field A and its conjugate A describe, respectively, the 
formation and destruction of pairs during the propagation. In particular, the definition of 
A makes use of correlated pairs of time-reversed single-particle states, in agreement with 
the original idea of Cooper |I3 ]. Generalizing slightly the development due to Gorkov [51 



we introduce such pairs by using an extended form of the time-reversed states, which we 
designate by ipT- Designating the time reversal operator by T, the usual time-reversed 
conjugate of the Dirac field operator ip is given by pE5|| 

^ T \x) = T^{x)7- X = B^ T (x) = 7 5^*(x) 

where 

x = (-£,£) , B = 7 5 C, 

and C is the charge conjugation operator. We define ipT as 

ifo(?) = r a ® 4) {r \x) = A^ T (x) , 

where A = r 2 <E> B and r 2 is the antisymmetric Pauli matrix, which operates here in the 
isospin space. Note that A = A T and A* = A^ . We then use the following ansatz for the 
effective single-particle Lagrangian 

JdtL eff = J d 4 x d 4 y fax) [ip-M + To//] 5(x - y)^{y) - ^(x)S(x, y)^(y) (4) 

+ \$(x)A(x,y)il)T(y) + -^ T {x)A{x,y)i(;{y) 

where 8{x — y) is a four-dimensional Dirac delta function and 

is the isospin matrix of chemical potentials, which will be used as Lagrange multipliers to 
fix the average number of protons and neutrons. 

The symmetries of the effective mean-field Lagrangian under transposition and Hermitian 
conjugation yield the following properties of the mean fields, 

A(x, y) = -A A T {x, y) A ] = -A (A{y, x)) T A ] and A(x, y) = -A A T {x, y) A\ (6) 
and 

= 7oS t (x,y)7o and A(x, y) = 7 A t (x, y)j ■ (7) 

The first of these symmetry conditions requires that the pair wave function be antisym- 
metric under exchange while the second guarantees real energy eigenvalues and probability 
conservation. 
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We can put the effective Lagrangian L e jf in a more symmetrical form by noting that 

d 4 x A ^(x) [ (ip - M + lof i)5(x -y)- y) ] ^(y) (8) 
= J d 4 xd 4 y i) T {x) [(ip + M — 7o/i)5(x - y) + T, T (x, y)] ^r(y) , 
where 

Z T (x,y) = Ai: T (x,y)A^ (9) 



The effective Lagrangian can then be rewritten in matrix form as 
1 



dtL e ff = - I d xd y (tp(x),ip T (x)) 



(10) 



(ip - M + ~iqh)8(x -y)- S(x, y) 
A(x,y) 



A(x,y) 



■>P(y) 



(ip + M -^ Q fj,)S(x-y) + E T (x,y) I \ tpr(y) 



which immediately yields the following coupled equations of motion for the fields ip and ipT, 
which we will call the Dirac-Gorkov equation, 



4 , (ifi-M + 'yon)5(x-y)-E(x,y) 



<?y 



A{x,y) 



4>(y) 



A{x,y) 



(ip + M -~fon)5(x-y) + Y,T{x,y) I \ ipr(y) 



= 0. (11) 



Defining the generalized baryon field operator as 



ip(x) 
ip T {x) 



we obtain a generalized baryon (quasi-particle) propagator 



S(x,y) 



G{x,y) F{x,y) 
F(x,y) G(x,y) 



m1) ) »M • 



(12) 



where, by (...), we mean the time-ordered expectation value in the interacting nuclear ground 
). We assume that the state 



state, (0 T( 



0) contains only nucleons interacting through 



the exchange of virtual mesons and contains no real mesons. 

We observe that G(x,y) is the usual baryon propagator while G(x,y) describes the 
propagation of baryons in time-reversed states. The off-diagonal terms of S(x, y) describe the 
propagation of correlated baryons and are the relativistic generalizations of the anomalous 
propagators defined by Gorkov |£3| . 

To derive the mean field equations, we first rewrite the interaction terms of the La- 
grangian density, Ci nt , as 



£-int{x\ 



x)r ja (x)<j)°*(x)i/j(x) 



(13) 



where the Greek letters «,/?,••• represent any indices necessary for the correct description of 
the meson propagation and coupling (Lorentz indices, isospin, etc.). The index j indicates 
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the mesons of the model: a, 00, it, and p, while their respective fields and meson-nucleon 
couplings are designated by the (j)J(x) and the Tj a (x). 



We then rewrite the meson fields 4>j in terms of their sources as 



<%(x) = Jd 4 yDf(x - y)i>{y)T iP {y)i>{y) , (14) 

where Dj l3 (x — y) is the Feynman propagator of meson j. Here, we have included the non- 
linear a-meson terms in the cx-meson propagator, so that they do not appear as corrections 
to the source term. Substituting in Eq. (p~3|) and inserting a factor of 1/2 for reasons of 
symmetry, we have 

JdtL int = -W j d A xd A y^{x)V ja {x)^{x)Df{x-y)^{y)V jP {y)^{y). (15) 
3 

Following Gorkov [44], we then obtain the mean field contribution of this interaction term 
by replacing each of the possible pairs of fermion fields by its vacuum expectation value, 

dt (L mt ) eff = E / d 4 xd 4 yDf(x - y) {*${x)T ja {x)ijj{x) @(y)r jf) (y)1>{y)) (16) 



+ 2iP(x)T ja (x) (iP(x)iP(y))T jp (y)iP(y) 

- ^{x)T ja (x) {^{xW{y)) T%(y)^ T (y) 

- V(x)Tj a (x) (^ T (x)^(y)) T jP (y)^(y) } . 

where (...) is again the time-ordered expectation value in the interacting nuclear-matter 
ground state. 

The first term in this expression is a Hartree one, the second a Fock exchange term 
while the last two, after using the definition of ipT to replace the transposed t/Vs, can be 
recognized as pairing terms. Comparing the mean field contributions to those of the effective 
quasi-particle Lagrangian, we can express the self-energy and pairing fields in terms of the 
two-fermion vacuum expectation values as 

E(x,y) = 6(x-y)Y,r ja (x) J d 4 z Df(x - z) (^(z)T^(z)^(z)) (17) 
j 

+ £ T ja (x)D^(x - y) U{x)~${y)) T jf3 (y) , 



and 

A(x,y) =Y,r ja (x)Df (x-y)(^(x)My))ATl(y)A\ (18) 
3 

while the equation for A(x — y) can be obtained from the equation for A{x — y) using the 
Hermiticity condition of Eq. (0). These expressions become self-consistency equations when 
we evaluate the expectation values by using their relationship to the generalized baryon 
propagator, Eq. (|i~2"D, which is itself a function of the mean fields. We find 
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y) = —i5(x — y)^2 ^ja( x ) J d^z D" /3 (x — z)Tr Y 'jp(z)G(z, z H 
j 



(19) 



and 



A(x,y) = T ja (x)Df(x-y)F(x,y)ATj p (y)A^ 



(20) 



The number of protons Z and neutrons N are the expectation values of the baryon number 
operators, N = i}){x)^q(\ ± r 3 )-^(x)/2, which we rewrite in terms of the generalized baryon 
propagator J/J, as 



Z 
N 



d 3 x ( ip(x)jo 



■ip(x) 



—i / d 3 xTr 



7o G{x, x 



(21) 



The Lagrange multipliers \i p and fi n , given in Eq. (H), are determined by requiring that these 
equations yield the desired values of Z and N. 

The Hamiltonian density operator is given by the T 00 component of the energy- 
momentum tensor, 



(22) 



Neglecting the retardation terms associated with the time derivatives of the meson fields 
in the ground state expectation value of Eq. (|22|), the energy density that results can be 
written as 



Six) =< H > =iTr 



(tf-d- M)G(x,x^ 



^g 3 a(x) 3 - ^g A a(xf 



(23) 



The total energy is obtained by integrating this density over space. 



III. PROPERTIES OF STATIC SOLUTIONS 

We will develop a static, ground-state solution to the self-consistency equations. We 
write the temporal Fourier transform of the full HFB propagator as 

v ' y ' ; ^ F(x, y; 00) G(x , y; u) ) ^ \ V a {x) j u - e a + irj V ayyh ayy 

+ (24) 
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The components U a ^ and V aj p are Dirac spinors corresponding to the normal and time- 
reversed components, respectively, of the positive-frequenc, e a , and negative-frequency, £p, 
solutions to the Dirac-Gorkov equation 

V A ^ y) (( £ - - y) + y)ho / V ^ / ' 

where we have introduced the single-particle hamiltonian, h(x,y), given by 

h(x, y) = {-ia ■ V + (5M)5{x -y)+ , y), (26) 

with 

hT(x,y) = Ah 7 (x^yjA^ and h(x,y) = (x,y). (27) 

After multiplying on the left by the matrix ^ ^° ^ V the Dirac-Gorkov equation, Eq. (P5D, 
can be written in Hamiltonian form as a Hermitian eigenequation, 

d 3 v ( {E + n)8{S-V)-h(S,jl) &(x,y) W C/(f) \ 

V A ( f ' ( £ - - v) + M*> £ ) / V t° v W J 

We thus conclude that the eigenvalues e are real. 

After multiplying the complex conjugate of Eq. Q2"5D on the left by the matrix 70 A ® t%, 
we can manipulate it into a form which is identical to the original equation, except for the 
sign of e. We may thus conclude that the solutions to the Dirac-Gorkov equation occur in 
pairs with real eigenvalues of opposite sign and eigenvectors of the form, 

e-e e--e .(^AV*(y)\ 

Using the above properties of the eigensolutions, we could rewrite the full HFB propagator 
of Eq. Q2~3|) as a sum over the positive eigenvalues e a alone. However, we will continue to 
use the form given in Eq. fl2"4|), distinguishing the positive and negative frequency solutions 
through explicit reference to one or the other. 

In the frequency representation, the self-consistency equations take the form 

E(f , y; u) = -i8{x -y)J2 ^a{^) J d 3 z Df{x - z; 0) J ^Tr [T l(3 (z)G(z, z; lo + )} (30) 

+*E/ ^r ja (f)^(f-y; W -^G(f,f;^)r j/3 (y), 

and 

A(x , M = 1 £ / ^T ja (x)Df(x -y;cu- u')F(x, y; u')ATj p (y)A^ , (31) 

where we have taken the vertices to be time independent and have assumed a time depen- 
dence of the form t — t' in all other quantities. We evaluate the equations in the static limit 
of the meson propagators, 
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D a/3 (x-y;uj) — ► D a \x - y; 0) = D a/3 (x-y). 
We then have for the self-consistency equations 



= -iS{x - y)Y^V ja (x) j d 3 zD" p (x - z) j — Tr \T jP (z)G(z, z; uj + ) 
j 

du 



(32) 



(33) 



+l Y^T ja {x)D^{x-y) J ^G(x,y;u) T 



jew) > 



and 



A(x,y) = iY,VU^Df{x-y) j j-F{^M A^0)A^ . 

3 

Finally, we evaluate the frequency integrals by closing the contour in the upper half-plane, 
yielding 



V{x,y) = 8($-y)Y i T ja (Z) j d 3 zDf(x-z) £ U,(z)T ]P (z)U,(z) 



(34) 



£t<0 



and 



A(5?,y) = -Er ja (x)< (a? -y) £ t/ 7 (*)K 7 (y) ^(y)^t 



(35) 



£t<0 



where the sum runs over the negative frequency solutions, e 7 < 0, of Eq. (|25|). 

The expressions for the number of protons and neutrons can be written similarly as 



Z 

N 



[ d 3 x J2 U^x)^-^-U 7 (x 



(36) 



E~,<0 



The expression for the energy can be reduced to 



E = Jd 3 x lj2U\(x)(e, + f,)U 7 {x)-^g 3 a(x) 3 -^ 

-i / rf 3 xrf 3 y £ (u\(xMx,y^U,(y) -U\(x)A(x,y)^ 



(37) 



Note that the last term is real due to the hermiticity of the Dirac-Gorkov equation. 

The summation over negative frequency solutions in Eqs. (|34]) through fl37D takes into 
account the occupation of all states, both those in the Fermi sea ( the surface of which is now 
diffuse, due to the pairing) and those in the Dirac sea. These expressions would thus require 
renormalization to yield finite results. Instead, we simply truncate the sums by excluding 
the solutions corresponding to quasi-particle states in the Dirac sea. This approximation 
has been found to give reasonable results in nuclear matter |29|. 
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Before reducing the self-consistency equations, together with the Dirac-Gorkov equation, 
Eq. (|25|), to a form in which they can be solved, we observe that, with B = 75C, each of the 




pair of wave vectors, I ^ )Z{ I and ( B ^ ^f] ) , possesses the time-reversed Dirac structure 



BV*{x) 

of the other. In Appendix A, we show that when these two states are equally occupied, then 
they satisfy the same Dirac-Gorkov equation with the same eigenvalue. We will assume this 
to be the case in the following. 



IV. REDUCING THE SELF-CONSISTENCY EQUATIONS 



We begin by analyzing the isospin structure of the self-energy and pairing fields under the 
assumption of pure proton-proton and neutron- neutron pairing. In this case, the solutions 
to the Dirac-Gorkov equation will be either purely proton particle-hole ones or neutron 
particle-hole ones of the form 



(U p \ 



\V P J 



and 



r 



V n 
\0 J 



where each of the elements in the column vectors are themselves 4-component Dirac spinors. 
Substitution of these into the self-consistency equations yields isospin structures of the mean 
fields of the form 



E{x,y) 



Z p (x,y) 

E n (x,y) 



and 



A{x,y) 



_ A p {x, y) 
A n {x,y) 



(38) 



where each of the elements of the matrices are themselves 4x4 matrices. The isospin de- 
pendent Dirac-Gorkov equation thus decouples into independent equations for neutrons and 
protons. In their Hamiltonian form, these are 



d 3 y 



(e + fi t )S(x -y)- h t (x, y) 



A]{x,y) 



u t (y) 

7oV t (y) 



A t {x,y) (e - Ht)5(x-y) +h t (x,y) y 

where we have written the respective Dirac hamiltonian operators as 

h t (x, y) = {-id ■ V + (3M)5{x - y) + /3£ t (x, y) t =p,n 







(39) 



(40) 



and the Lagrange multipliers as \i v = fi + 5fi and \i n = /i — SfM. We have also made use of the 
fact that T^tr{x, y) = T> t {x, y), so that h t r{x, y) = h t {x, y), since we have assumed invariance 
under time inversion of the Dirac structure. 

We want to obtain the coupled equations for the case of an axially-deformed nucleus. 
We take the z-axis to be the symmetry axis and use cylindrical coordinates {x = r^costp, 
y = r± sin (p , and z) . Although the total angular momentum j is no longer a good quantum 
number, its projection along the symmetry axis Q, as well as the parity it, (and, of course, 
the isospin projection, which we also denote by t) are still good quantum numbers. We write 
each of the 4-component spinors of the wave-function as 
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U tl {x) 



uf tJ (x) 
iug tl {x) 



and 



vft-y(x) 
ivg tj (x) 



uf t j(r ± ,z) e l(n - +1 / 2 ^ 
iugty(r±,z) e^- 1 ' 2 ^ 
\iug^(r ± ,z) e l ^' +l l 2 ^ ) 



I vf^r^z)^^- 1 ' 2 ^ \ 
vf t -(r ± ,z) e l{n -< +1 ' 2 ^ 
ivg£ f (r ± ,z) e i{ ^~ 1 / 2 ^ 

\ivg^(r ±: z) e i{n ~< +1 ' 2 ^ J 



(41) 



(42) 



in which the functions uf^, ug^, vf^ : and vg^ are real. For each solution with positive value 
of the angular momentum projection, f2 7 , we have a time- reversed solution with the same 
energy but a negative value of the angular momentum projection, -f2 7 , given by BU^(x) 
and 7o-Bl / t *(a 7 ). The densities that enter the Hartree terms of the self-energy can then be 
written as 



Ps(r±,z) 
Pb(t±,z) 
Ps(r±,z) 
Pc(r±,z) 



2 £ UtrloUpr, 

w t7 <0,t,C 7 >0 

2 £ Ulu tl , 

u>t~/<0,t,Qj>0 

2 £ 2m t UlU tJ , 

w t7 <0,t,C 7 >0 

2 E K + V2) c/ t f 7 c/t 7 , 

w t7 <o,t,n 7 >o 



(43) 



in which the sum over states of different parity is implicit in the sum over states of different 
energy. The local Hartree contribution to the self-energy may be written in terms of these 
densities as 



f3Y, H (x,y) = S(x - y) j d 3 z (- (3 g 2 a d a (x - z) p s (z) + gld° w (x - z) p B (z) 



(44) 



+ 



T 3 d°(x - z)p 3 (z) +e 2 



d Jx-z)p c (z) , 



where the meson propagators have been reduced to 

1 exp (— rrij \x — y\) 



dPAx-z) = 



An 



\x - y\ 



(45) 



with the exception of the a-meson propagator, which is still assumed to contain the nonlinear 
terms. It is customary to write the Hartree contribution to the self-energy in terms of the 
mean fields associated with each of the mesons, as 



^ H (x) = -(3g a a(x) + g u u\x) + | r 3 p 00 (x ) + e A °(x) , 



(46) 



with 
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uj°(x) = g w j d 3 zdl(x- z)p B (z) 



p 00 (x) = ^Jd 3 zd° p (x-z)p 3 (z), 

A°(x) = e f d 3 z d°Jx - z)p c (z) , (47) 



a{x) = g a J d zd a (x- z)p s (z) 

= J d 3 z d° a (x - z) [g a p s {z) - g 3 a(x) 2 - g 4 a(x) 

In the last expression, we have written the a mean field in terms of its free propagator, 
as given in Eq. (ffH), and included the nonlinear terms explicitly. The mean meson fields 
possess the same axial symmetry as the densities. 

The Fock exchange term of the self-energy has the form 

e f (x, y) = -J2 ?U?) DaP ^ - V) (48) 

3 

x £ Xtxl (Ut,(x)uUy) + BU^(x)UT(&)Bl) l0 T^(y) , 

e t7 <0,i,f2 7 >0 

where we have written the isospin dependence of the intermediate states in terms of the 
isospinor Xt- As has been discussed in the literature (cf. Ref. or Ref. f42||), the most 
important effects of the Fock terms due to the exchange of the short-range a, u, and p 
mesons can be taken into account by using adjusted Hartree terms. This can be better 
understood by looking at the zero-range limit of the meson propagators, for which the Fock 
terms reduce to linear combinations of the Hartree ones. Here, we retain the small but finite 
range of the meson propagators, but approximate the self-energy as the Hartree contribution 
alone and use the parameters appropriate for this purpose. 

The pairing field arises from an exchange term, similar in form to the Fock term of the 
self-energy. We can write it as 

A(x, y) = -J2 r ja (x-)Df(x - y) (49) 
j 

x E XtXi (u tl {x)V t \{y) + BU^{x)V^{y)B^) lo AY%{y)A\ 

et 7 <0,t,f2 7 >0 

where xp,n = Xn, P are the time-reversed isospinors. 

We neglect the Coulomb contribution to the pairing field here and approximate the 
contributions of the other mesons using the zero-range (infinite-mass) limit of the meson 
propagators. We also neglect the contribution of the non-linear cx-meson terms. As we have 
argued in the introduction, due to the short range of the exchanged mesons, a zero-range 
approximation to their propagators does not make a significant difference in the results. It 
also greatly simplifies the numerical calculations. The Hamiltonian form of the pairing field 
is, in this case, 

Al(x,y)= lo M^yho (50) 
= 5{x - y)c pair (^lo Kt{x) 7o - + ^ ) «t(x) 7m7o) , 



14 



where the anomalous density n t {x) is 

n t {x) = (Ut 7 (x)V^(x) + BU; 7 (x)V t T(x)Bi) 70 

£t 7 <0,f2 7 >0 



t = p,n. 



(51) 



However, the pairing field, in the given form, does not necessarily satisfy the antisymmetry 
condition of Eq. (||). Applying the antisymmetry condition to the pairing field, we can 
reduce it to a condition on the anomalous density, 



K t (x) = BkJ(x)B^ . 
To ensure that this is satisfied, we take the anomalous density to be 

K t(x) -> - (n t (x)+BnT(x)B^ 



t — p,n 



2 

= \ E (u t7 (x)V^(x) Jo + M/ t ;(£)^(£)5t 7o 

A £ t7 <0,C 7 >0 

+ l0 V tl (x)U} 1 (x) + l0 BV t ;(x)U^x)B^) . 
We then write the components of the pairing field as 



&l(x) 



( 5it(r±,z) iS 3t (r±,z) i 5 4t (r±, z)e~ iLp \ 

6 u {r_L,z) i5 u (r±,z)e llp -i5 3t (r ± ,z) 

-i$3t(r±,z) -i5 u {r L ,z)e- lLp 5 2t (r±,z) 

\-i$4t(r±,z)e i ' p i5 3t (r±,z) S 2t (r ±J z) J 



where the four real functions 5j t are 



Su(r±, z) = (c s - c v )K U (r±, z) - 3c v K 2t (r ± , z), 

S 2 t{r±, z) = (c s - c v )K 2t (r±, z) - 3c v K lt (r ± , z), 

$3t(r±, z) = (c s + 2c v )K 3t (rj_, z), 

S*t(r±, z) = (c s + 2c v )k u (t±, z), 



with 



■^pair n 

mi 



and 



pair 



' gj , {qM 

m 2 , m\ 



The four real components of the anomalous density are 
Ku(r±, z)= E (ufttfVf^ + uf t ~vf; 



K2t{r±,z) 
K>3t{r±,z) 
K4t(r±,z) 



't-yji 

e t7 <0,f2 7 >0 

£ t7 <o,n 7 >o 

et 7 <0,f7 7 >0 



£t 7 <0,O 7 >0 



(52) 



(53) 



(54) 



(55) 



(56) 



(57) 
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The Dirac structure of the pairing field is very similar to that of the 1 So pairing field in 
symmetric nuclear matter, where 

A nm (k) = 7 A nm t (fc)7o = A s (k) - 7 A (k) - ry 7 ■ k A T (k) . (58) 

The upper and lower diagonal components o fthe pairing field, 5\ t and §2t can De directly 
associated with the linear combinations As ± Ao of the nuclear matter components, while 
the remaining components, 8%t and Su are more loosely related to the contributions of the 
nuclear matter tensor term, At- 

An overall constant c pa i r has been introduced in the expression for the pairing field to 
compensate for deficiencies of the interaction parameters and of the numerical calculation. 
The necessity for such a constant is apparent from studies of pairing in nuclear matter. Non- 
relativistic [{46-48|| and relativistic [M] calculations have verified that 1 Sq pairing in nuclear 



matter is dominated by the two-nucleon 1 So virtual state. Pairing in nuclear matter is weaker 
the further the x So virtual state is from the real axis in the complex-momentum plane. The 
location of the virtual state depends on the strength and form of the two-nucleon interaction 
and on the space of states used in the calculation. In Ref . , various sets of interaction 



parameters, even zero-range ones, were shown to furnish mutually consistent physical values 
for the pairing gap function, when they were supplemented with a large momentum cutoff 
adjusted so as to place the two-nucleon virtual state at its physical location. We expect 
a condition similar to that in nuclear matter to apply here. However, as it is extremely 
difficult to fix the position of the two-nucleon virtual state within the harmonic oscillator 
basis that we use, we instead multiply the pairing field by an overall constant that we expect 
to be able to fix independently of the charge and mass of the systems under consideration. 
We emphasize that this is not a weakness of our calculations alone, but of any Hartree- 
(Fock)-Bogoliubov calculation using a limited space of states and an effective interaction, 
even those using a finite-range one. The pairing field obtained in such a calculation will 
depend on both the interaction and the space of states used and will usually require that 
one or the other of these be adjusted on order to obtain reasonable results. Here, we find it 
more convenient to introduce an arbitrary constant in the interaction rather than arbitrarily 
limit the space of states we use. 

With the above simplifications in the self-energy and pairing fields, the Dirac-Gorkov 
equations for neutrons and protons reduce to local differential equations. Their Hamiltonian 
form is 



e + Hj-h t (x) A\(x) 
A t (x) £ - /i t + h t (x) 




0, (59) 



with 



where 



h t (x) = -ia ■ V + f3M*{x) + V t (x) t = p,n, (60) 



M*(x) = M-g a a{x), (61) 
V t (x) = g u uo\x) + 1 2m t p 00 (x) + e (1/2 + m t ) A°(x) , 
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and Aj(x) is given in Eq. (|54]) . 

The total energy can now be written in terms of the mean fields as 



E = J d?x ( £ |t/ 7 (f)| 2 ( £7 + /i)-i(73a(x) 3 -i 

i & p 00 (f ) P3 (x) - \e A\x) p c {x) + ~ £ Tr[Al (f ) ^ 



^ 4 ^(^) 4 + -g a a(x) p s (x) (62) 



In the expression above, we have also subtracted the harmonic oscillator estimate to the 
center-of-mass motion, 

E cm = hu = MeV, (63) 

in order to obtain an expression for the total internal energy of the nucleus. 



V. NUMERICAL SOLUTION OF THE DHB EQUATION 

We solve the Dirac-Gorkov and the Klein-Gordon equations by expanding the fields as 
well as the wavefunctions in complete sets of eigenfunctions of harmonic oscillator potentials. 
In actual calculations, the expansion is truncated at a finite number of major shells, with the 
quantum number of the last included shell denoted by Np in the case of the fermions and by 
Nb for the bosons. The maximum values are selected so as to assure the physical significance 
of the results obtained. The same procedure has been used by many researchers, among 



them, by Vautherin |49| in the nonrelativistic Hartree-Fock approximation, by Ghambir et 
al. [IT] in the relativistic mean field + BCS approach and by Lalazissis et al. [f40|-[42|l in the 
RHB approach. 

The spinors of the Dirac-Gorkov equation are expanded in terms of the eigenfunctions 
of an axially-deformed harmonic-oscillator potential, 

V osc (r ± , z) = X -Mulz 2 + ^Mulr 2 . (64) 

The oscillator constants are taken as 



with volume conservation relating the two constants to that of a spherically-symmetric 
potential b\b z = b%. 

The eigenfunctions of the deformed harmonic oscillator can be written explicitly as, 

gj mi ip 

$a(r) = C'K) V>n.(z) -^=Xm a Xra t (66) 

V Z 7T 

where a denotes the complete set of quantum numbers (n r , mi, n z , m s , and m t ) and 
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N mi / r \ 2 

<J(r±) = -£-y/2rr /2 L%(Tj)e-*'' a with 77 = J (67) 

^W = ^n,(Oe^ 2/2 with 



In Eq. (|67|) , L™*{rj) and H nz {^) are Hermite and associate Laguerre polynomials fl50l , with 



the normalization constants, N™J and N Hz , given in Ref. |T7| . In these equations, n r and 
n z are the number of nodes in the r and z directions, and mi and m s are the projections 
of angular momentum and spin on the z axis. The third component of the total angular 
momentum f2 7 and the parity 71 are then defined as, 

fi 7 = m { + m s 71 = (-l)^+ m < . (68) 

We expand the Pauli components of the Dirac spinors, uft 1 (r±, z), ug t -y(r±, z), vf tl (r±, z), 
and vgt~f(r±, z), in terms of the oscillator eigenfunctions. Inserting these expansions into the 
Dirac-Gorkov equation Eq. (|59|), we can reduce the equation to the diagonalization problem 
of a symmetric matrix and calculate the Hartree densities of Eq. (f43l) and components of the 
anomalous density of Eq. (|57|) The fields of the massive mesons are expanded in a manner 
similar to the fermion expansion, with the same deformation parameter f3 but a smaller 
oscillator length of b B = b /\^2. The Coulomb field is calculated directly in configuration 



space. In short, the method used is a direct generalization of that described in Ref. |L7 
where more details may be found. 



VI. NUMERICAL RESULTS 

The parameters required to perform numerical calculations are the nucleon and meson 
masses, the meson-nucleon coupling constants and the factor c pa i r that multiplies the pairing 
interaction. The calculations that we present were performed using the masses and coupling 
constants of the NL3 potential ||51|| . We performed calculations for several values of the 



factor Cp a i r . These permit us to study the extent to which physical observables depend on 
the pairing interaction and to choose the value of the parameter c pa i r that best fits the 
observables. 

As the mesons fields and the nucleon wavefunctions are expanded in a deformed basis of 
harmonic oscillator states, we must also specify the number of major oscillator shells to be 
used in the expansions for fermions, Np, and bosons, Nb, as well as the basis deformation 
parameter, (3q. Here, we work, for the most part, in a basis of 12 major oscillator shells for 
fermions and 24 for bosons. We use the standard expression for the oscillator frequencies, 
Tiujq = 41A -1 / 3 MeV, and spherical bases, with (3q = 0, in all calculations. Although the 
number of major oscillator shells is not as large as might be desired, it does seem to be large 
enough to obtain reasonable values for the observables studied. 

In order to analyze the characteristics of of the single-particle levels, we define several 
related average values. First, in deference to the standard nonrelativistic notation, we define 
the occupation probability of each of the two states of the level with frequency e tl < as 

v% = f d 3 x\U tl (x)\ 2 , (69) 
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so that we also have 



u 



1-v 



try 



d 3 x |V^ 7 (£)| 2 . 



(70) 



(71) 



in which we take advantage of the normalization of the state vector to avoid normalizing 
the result. This is not possible for the pairing term, for which we define the gap parameter 
of a single-particle level as 



We define the energy of a single-particle level as 

E tJ = [ d 3 x (U^(x) h t U tl {x) + V t \{x) h t V tJ (x) 



A 



'- 



d 3 x U^(x) Al<y Q V tl (x) / u tl v tT 



(72) 



The definition of the gap parameter is fragile and is subject to numerical error when t> <7 — > 1 
or Ut-y —* 1. We define an average gap parameter for the neutrons and protons of a nucleus 
as 



< A t >= A 



*7 U t-y 



£-v<0 



£-,<0 



y *7 



(73) 



We note that, unlike in the BCS approximation, no relation exists here between the oc- 
cupation probabilities, i>f and w 2 7 , the energy E tl , and the gap parameter A t7 . However, 
these averages still furnish a good description of the most important characteristics of the 
single-particle levels. 

In the following, we first analyze pairing in the isotopes of tin and then discuss similar 
results for the isotopes of nickel and calcium, all of which are spherical. The isotopes of tin 
and nickel have already been the object of two thorough studies ||i0|jn^ . We then analyze 
deformation and pairing of the N=28 isotonic chain, which has been the object of two 
recent studies P2|,p8|,|59| . Finally, we turn our attention to nuclei in the region of the Z w 40 



subshell closure, the isotopes of Kr and Sr, in particular. These have been the object of 
many previous studies. p0|-E2,El,EH 



A. Spherical Nuclei 

Pairing in the Sn and Ni isotopes has been studied extensively, due to the simplification 



provided by their spherical symmetry. In Ref . |52| , Dobaczewski et al. describe a very com- 
plete study of pairing in the Sn isotopes, based on a nonrelativistic Hartree-Fock-Bogoliubov 
approximation, in which they examine the differences in pairing due to the use of various 
effective interactions and due to coupling to the particle continuum, as well as the effects of 
these on experimental observables. In Ref. j|D), Lalazissis et al. compare RHB calculations 
of both the odd and even isotopes in the Sn and Ni chains with the experimental data, 
obtaining quite good agreement. In another study Meng compared RHB calculations 



using the finite-range Gogny D1S and a zero-range interaction and found good agreement 
between the two and with the experimental data. 
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We begin our study of spherical nuclei with the tin isotopes. We performed calculations 
of the ground states of the even isotopes from 100 Sn to 176 Sn, that is, from the closed neutron 
shell at N=50 to the closed shell at iV=126. We present calculations for three values of the 
parameter c pa i r . The objective of this study was to identify the observables sensitive to the 
parameter c pa i r (and, thus, to the pairing) and to adjust the parameter accordingly 

We present in Fig. [I] the two-neutron separation energy of the even Sn isotopes in the 
mass range from v4=100 to A=176, calculated for the values c pa i r =0A5, 0.47 and 0.50. We 
compare our calculations with the two-neutron separation energies obtained using the ground 
state masses tables of Moller-Nix and of Audi-Wapstra J57| . The Moller-Nix values for 
ground state masses were calculated using an extended finite-range droplet model with 
parameters adjusted to the experimental ground state masses. The Audi-Wasptra values for 
the ground state masses are essentially the experimental ones, with some extension to proton- 
and neutron-rich nuclei based on systematics. We verify that all three calculations follow 
the trend of the Moller-Nix and Audi-Wapstra values. The two-neutron separation energy 
is seen to be almost independent of the strength of the pairing field, with the calculation 
using c pa j r =0.50 providing only slightly better agreement with the data than the others. 
All the calculations underestimate S2n m the region of v4=100. This discrepancy might 
be attributable to a deficiency in the isospin dependence of the NL3 parameter set. It 
could also be due to effects that are not included in the calculation, such as neutron-proton 
correlations, which are suspected of being of importance in N « Z nuclei. Our two-neutron 



separation energies are in good agreement with the RHB ones of Ref. ^0|. However, such 
consistency is to be expected, given the relative insensitivity of the two-neutron separation 
energy to the pairing, since both calculations use the NL3 parameter set. 

The calculated two-neutron separation energies remain almost constant from the shell 
closure at A=132 to about A=160, in contrast to the Moller-Nix values, which decrease 



slowly. The SIII and Sk5 Skyrme interaction calculations of Ref. [[5^] also furnish a relatively 
constant separation energy above the v4=132 shell closure, which extends to even higher 
values of the mass, although most of their calculations are in agreement with the Moller-Nix 
systematics. The nucleus 176 Sn is unbound in the calculation with c pa j r =0.50, as it is in all 



of the calculations of Ref. fl5"2"| . In our calculations, this is due to the fact that 174 Sn is more 
tightly bound for c pa j r =0.50 than for the other values of the parameter. The binding energy 
of the magic nucleus 176 Sn is the same for all values of the parameter c pa i r - 

In Fig. ^| we show the average value of the neutron gap parameter of the even Sn isotopes 
as a function of the mass number A, for the same three values of c pa i r . The gap parameter 
possesses a clear dependence on the value of c pair . We also show in the figure the Moller- 
Nix and experimental Audi-Wapstra values for the standard estimate of the neutron gap 
parameter as the difference between the binding energies of an even-even nucleus and its 
odd mass neighbors, 

< A n (Z, A) >= B(Z, A) - \{B(Z, A - 1) + B(Z, A + 1)) 

where B(Z, A) is the binding energy. The calculation with c pair =0.50 shows reasonable 
agreement with the Moller-Nix values in the region of the N=82, A=132, shell closure, but 
tends to underestimate the Moller-Nix values well below the shell closure and overestimate 
these values well above the shell closure. 
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The neutron shell closures at iV=50,82, and 126 are clearly visible in Fig. |2|. At each 
of values of the neutron number, calculated neutron gap parameter goes to zero. We note 
that the experimental gap parameter has maximum, rather than a minimum, at each shell 
closure, reflecting the local maximum of the binding energy that occurs there. Several 
subshell closures are also visible in the calculations, at values of the mass at which the gap 
parameter reaches a non-zero minimum. These occur at N=58, A— 108, between the l<?7/2 
and 2<i5/2 levels, at N=Q4, v4=114, between the 2<i5/2 level and the remaining levels of the 
4hu shell, and at iV=112, A=1Q2, between the 5hui shell and the 1H3/2 level. Subshell 
closures are more visible when the pairing interaction is weaker and, thus, more sensitive to 
the energy differences between the levels. 

Average values of the neutron gap parameter of the even Sn isotopes obtained in RHB 
and nonrelativistic calculations were presented in Refs. ffl] and |p2 |, respectively. The 



calculations using the Gogny interaction, both the RHB and the nonrelativistic ones, furnish 
a large value of the gap parameter - above 2 MeV - and show no subshell structure. Were 
we to increase the strength of the pairing interaction so as to obtain similar values for the 
pairing gap, our calculations would also show no shell substructure in the average pairing 
gap. The nonrelativistic calculations of Ref. ||5"2" |, using Skyrme pairing interactions produced 
results in closer agreement with ours. These calculations also show structure due to subshell 
closures, although the structure is different from that seen here, possibly due to differences 
in the spin-orbit splitting of the levels, due to the different mean fields. 

In Fig. 0, we show the pairing energy of the even Sn isotopes as a function of the mass 
number, obtained using the same three values of c pa i T . The pairing energy is defined as 



where K t and A] are given in Eqs. (|53|) and (|54"|) , respectively. We observe that the pairing 
energy E pr displays a dependence on the parameter c pair and a shell structure similar to 
those seen in Fig. [| for the average gap parameter. However, E pr is observable only through 
its effects on the binding energy, which diminishes its utility as a means of determining the 
parameter c pa i r . 

In Fig. f| we present the difference between the calculated value of the binding energy of 
the even Sn isoptopes and the Moller-Nix value, as a function of the mass, using the same 
three values of c pair used before. We also present the difference between the calculated value 
for c pa j r =0.50 and the Audi-Wapstra values, where the latter exist. In the mass range of 
the valley of stability, the difference between the two binding energies can be minimized, to 
a certain point, through an appropriate choice of c pair . For the tin isotopes, a value of c pair 
close to 0.50 seems to yield the best average agreement. A value much larger than this would 
destroy the good agreement obtained for the gap parameter of Fig. |[ However, a larger 
value of c pa i r would smooth out the two dips in the binding energy difference that roughly 
follow the mass dependence of both the average gap parameter and the pairing energy. This 
would yield pairing gaps in better agreement with the Gogny D1S ones of Ref. We 
then might expect to be able to reduce the remaining discrepancy by adjusting the mean 
field parameters. 

In Fig. |], the calculated value of the binding energy is larger than the experimen- 
tal/systematic values in the N « Z region and reaches a difference of about 4 MeV for 
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the magic nucleus 100 Sn, a result which is independent of the value of c pa i r . The sharp 
increase of the binding-energy difference for iV « Z is consistent with the discrepancies be- 
tween the calculated and the Moller-Nix values of the two-neutron separation energy, shown 
in Fig. |l|. We can attribute the discrepancies to deficiencies in the isospin dependence of 
the NL3 parameters. If we were to include neutron-proton pairing in the calculations, the 
difference would only increase as the pairing would bind even more the already overbound 
N fa 50 nuclei. The large differences in binding energy on the neutron-heavy side of the 
curve could be due as much to deficiencies in the Moller-Nix systematics as in the NL3 
parameters. 

In Fig. |5| we present the deviation from the systematic value, r^N 1 ^ 3 , of the root mean 
square neutron radius of the even isotopes of Sn, as a function of the mass. The deviation 
displays a clear minimum at the principal shell closure at N=82. Although the deviation 
also tends to decrease near the shell closures at iV=50 and 7V=126, the effect of these more 
extreme shell closures is much smaller than that of the N=82 one, within the stability 
valley. The deviation from the systematic value also possesses a sharp maximum at iV=112, 
between the 5hu shell and the 1*13/2 level. We recall that the average gap parameter and 
the pairing energy display a minimum at the same subshell closure. Other subshell closures 
also appear in the deviation of the radius, but only as changes in its slope. The deviation 
from the systematic value of the root mean square mass radius of the even Sn isotopes 
shows the same structure, but to a lesser degree, and is not shown here. The deviation from 
the systematic of the root mean square proton radius, which is also not shown, increases 
monotonically with the mass number. 

We find that the deviations of the root mean square radii are almost independent of 
the parameter c pa i r and, thus, of the pairing interaction. Slightly larger variations have 
been seen in a comparison between RHB and Hartee+BCS calculations El. Still, with the 



exception of a few very special nuclei, such as 11 Li, the root mean square radii of a nucleus 
would seem to be determined almost exclusively by its mean field. 

As the second example of a spherical isotopic chain, we studied the even isotopes of 
nickel, performing calculations from 48 Ni, at the iV=20 shell closure, to 100 Ni, two neutrons 
beyond the subshell closure at iV=70. We again performed calculations for three values 
of the parameter c pa i r , their values in this case being c pa j r =0.49, 0.50, and 0.52. Here we 
wanted to verify the generality of the observations made in the case of the Sn isotopes. 

The calculations of the two-neutron separation energy for the Ni isotopes describe the 
tendency of the data fairly well, reproducing the discontinuities in separation energy that 
occur at the shell closures at N=28 and iV=50. The calculated values of the (two-neutron) 
separation energy are almost independent of the parameter c pair and underestimate the 
mass dependence of the separation energy at very low and very high values of the neutron 
excess, as in the case of Sn. Not surprisingly, our values are in good agreement with the 



RHB ones of Ref. [|0|, which also use the NL3 parameter set. They also agree well with 
the results of Ref. |3Sj| , obtained with the NLSH parameter set |53| , using Gogny D1S and 
density-dependent zero-range pairing interactions. 

In Fig. H we present the mean value of the neutron gap parameter of the even Ni isotopes 
as a function of the mass, for the same three values of c pa i r - Again, we find a clear depen- 
dence of the gap parameter on the value of c pa i r - The calculation with c pa j r =0.52 shows the 
best agreement with the Moller-Nix and Audi-Wapstra values in the region above the shell 
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closure at iV=28, A=56. All of the calculations are in disagreement with the experimen- 
tal/systematic values in the region of the magic number iV=28. This again suggests that 
shotr-range neutron-proton correlations may be important when N pa Z. 

In Fig. |6], we note that the neutron gap parameter is zero at the neutron shell closures at 
iV=20, 28, and 50. The subshell closures are more marked here than in the case of tin. For 
Ni, these appear at N=40, A=68, between the 3hu shell and the lg9/2 level and at N=70, 
A=98, between the 4hu shell and the lill/2 level. As occurred in the calculations of the 
isotopes of Sn, the value of the average gap parameter at the subshell closures increases with 



Cpair • 



Values of the average neutron gap parameter of the Ni isotopes between iV=28 and 
iV=50, obtained using the Gogny interaction in a RHB calculation, were presented in Ref. 
iOfl . As in the case of Sn, the RHB calculation furnishes values of the gap parameter greater 



than 2.5 MeV, which are quite large. The calculation presents a subshell structure similar 
to that with c pair =0.52 in Fig. || but almost 50% larger in magnitude. 

As in the case of the Sn isotopes, the pairing energy of the even Ni isotopes reflects 
the shell structure and the dependence on c pa i r observed in the average gap parameter, but 
furnishes no new information. The root mean square radii reflect the same shell structure, 
although with less clarity, and are almost independent of the value of c pa i r . 

We also examined the Ca isotope chain, performing calculations of the even isotopes of 
Ca from N=8, A=28, to iV=50, A=70, for two values of the parameter c pa i r . Comparisons 
between the two calculations and the Audi-Wapstra and M611er-Nix systematics resulted 
in conclusions similar to those obtained for the tin and nickel isotopes. In particular, the 
two-neutron separation energies were well described in the two calculations, which show 
almost no dependence on the parameter c pair . A value of c pa j r =0.55 best agrees with the 
gap parameters obtained from the Audi-Wapstra systematics, which, in this case, lie about 
20% below the Moller-Nix values. 

In summary, we found our calculations of the root mean square radii and the two-neutron 
separation energies of spherical nuclei to be relatively independent of the strength of the 
pairing interaction. As the experimental data for these quantities are well fit by the NL3 
parameter set used here, it is not surprising that our calculations describe them well. We 
observed that the pairing energy is only observable through its effects on the binding energy 
and that our calculations are not as successful at describing the latter as they are with the 
radii and the separation energies. A reasonable fit to the binding energies will depend on 
the adequate choice of the parameter c pa i r as well as further adjustments in other interaction 
parameters. We found that the binding energy and the mean value of the gap parameter, 
to a certain point, are the observables that can be used as guides to a study of pairing in 
spherical nuclei. 



B. Several deformed light nuclei 

In this section, we analyze pairing and deformation in the N = 28 isotonic chain. The 
measured quadrupole deformations for these nuclei characterize 48 Ca as spherical but 46 Ar 
and 44 S as deformed, implying suppression of the N = 28 magic number. The deformation 
can be explained, in this case, by the close proximity of the 1/7/2 level to the 2p levels in 
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the 4hu shell. In the nucleus 48 Ca, the spherical proton configuration constrains the neutron 
configuration to also remain spherical. In the 46 Ar and 44 S nuclei, the partially filled proton 
shell perturbs the neutrons sufficiently for them to prefer a deformed configuration in which 
the 1/7/2 and 2p orbitals are partially occupied. 

Recent theoretical studies describe the experimental data for these nuclei reasonably 
well with nonrelativistic and relativistic mean field + BCS |58| , |59| calculations and with 
RHB calculations using the Gogny potential J4"2" |. In Refs. [5S| and [ |5"9[] , Werner et al. use 



both a Skyrme-Hartree-Fock approach and a relativistic mean field model, with the NLSH 
parameter set, to study various istopic chains in the N « 28 region. They include pairing 
as a small (75 keV) constant pairing gap in a BCS formalism. They found the isotones 
42 Si, 44 S and 46 Ar to show evidence of shape coexistence in both approaches, although the 
deformation at the energy minimum and its energy difference with the excited minimum 
varied between the two calculations. Lalazissis et al. used the RHB approach to perform 
a similar study of the iV « 28 isotopes in Ref. ||42|| . Using the NL3 parameter set, they 
studied the even isotopes of Mg through Cr. They observed prolate-oblate staggering in the 
iV = 28 isotones below 48 Ca, similar to that seen in the RMF calculations of Werner et al. 
and similar signs of shape coexistence: the binding energies versus quadrupole deformation 
of the isotones display two minima or a very flat single minimum. 

We present here the average values of the gap parameter of the even-even S isotopes, in 
Fig. [7| as a function of the mass, for two values of the parameter c pair . We observe that both 
calculations fall well below the Moller-Nix and Audi-Wapstra values at lower values of the 
mass. However, the isotopes that are of interest at the moment are in the region of large 
neutron excess, A>40, where the average values of the calculated gap parameter and the 
experimental/systematic values are in reasonable agreement. In particular, we see in Fig. [7], 
that the Moller-Nix and Audi-Wapstra values of the gap in 44 S fall close to the calculation 
with Cp a i r — €.55. The average value of the gap parameters of the other isotopes that we will 
examine show the same reasonable agreement with the experimental/systematic values in 
the mass range of interest. 

To study in detail the dependence of the binding energy on the quadrupole deformation, 
we perform calculations at several fixed values of the deformation. To do this, we include in 
the Lagrangian a term that is quadratic in the difference between the quadrupole moment 
of the nucleus and the desired value of the moment. The solution of the equation of motion 
including this additional potential term tends to take a value which minimizes its contribu- 
tion, thereby yielding a quadrupole moment close to the desired value. The contribution of 
the constraint term to the energy is subtracted to obtain the binding energy of the system 
at the value of the quadrupole moment obtained in the calculation. The consistency of the 
method can be verified by calculating the ground state energy of a spherical nucleus as a 
function of the quadrupole deformation (3. For the case of 48 Ca, we find that the ground 
state energy reaches its minimum at (3 = 0, where the nucleus is spherical, just as we would 
expect. We also find a small difference of about 0.5 MeV between the calculations for the 
two values of c pa i r . This difference is consistent with the small difference in the pairing 
energy and the small, but nonzero, average pairing gap found in the calculations. 

In Fig. |8|, we present the dependence on the deformation (3 of the binding energy of the 
nucleus 44 S, for the two values of c pa i r . (The energy Ej, displayed in the figure actually differs 
by a sign from the binding energy, Eb(Z,A) = —B(Z, A).) Here the variation of c pa i r from 
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0.50 to 0.55 makes a difference of almost 3 MeV in the binding energy, in contrast with the 
small difference of about 0.5 MeV found for 48 Ca. The energy curves for 44 S possess two 
well defined minima, with the oblate minimum about 1 MeV above the prolate one in the 
case of c pair =0.50 and about 0.5 MeV above the prolate minimum when c pair =0.55. The 
increase in the parameter c pa i r tends to flatten the energy curve, thereby diminishing the 
differences between its peaks and valleys. The increase in c pa i T also reduces the value of the 
deformation at the minima. The deformation at the oblate minimum decreases in magnitude 
by almost 20% with the increase of c pa i r from 0.50 to 0.55 , while the deformation at the 
prolate minimum is reduced by about 10% with the same increase in c pair . 

The experimental value of the quadrupole deformation of 44 S is (3— 0.258(36) pO| , |6T| , 
which is in good agreement with the calculated value, /3=0.28, using c pa j r =0.55. The rela- 
tivistic calculations of Werner et al. and Lalazissis et al. furnish a value of the deformation 
similar to this one. The RHB calculation of Ref. |42| yields an oblate minimum about 200 



keV above the prolate one, which can be compared to our value of about 0.5 MeV between 



the two minima and a value of about 0.8 MeV in the RMF+BCS calculation of Refs. 58 



and [ESI] . Nonrelativistic HF calculations furnish a value of (3 that is about half of the 



experimental one [|58,|59]. 



In Fig. we present the dependence on the deformation (3 of the binding energy of the 
nucleus 46 Ar, for two values of the parameter c pa i r . The difference in the binding energy due 
to the variation of c pa ; r is about 1 MeV. The minima of the curves are much less defined here 
than they are in the case of 44 S. In fact, the calculation using c pa j r =0.50 appears to have only 
an oblate minimum. The calculation with c pa j r =0.55 possesses two minima, with the oblate 
one approximately 200 keV below the prolate one. Due to the lack of structure in the curves, 
it is difficult to determine the reduction of the deformation at the minima due to the increase 
in the parameter c pair . The experimental value of the quadrupole deformation of 46 Ar is \(3\ = 
0.176(17) fl6~2"l , with the sign of the deformation undetermined experimentally. Assuming the 
deformation to be oblate, the calculated value of /3=-0.18 is in excellent agreement with the 
data. A nonrelativistic HF calculation furnishes results in reasonable agreement with ours, 
while its companion RMF+BCS calculation predicts the oblate minimum of the 46 Ar energy 
curve to lie slightly above a spherical ground state p8| , |59| . The RHB calculation using the 
Gogny interaction yields a flat energy curve with a barely discernible minimum at f3 ^-0.15 



We have performed calculations of the binding energy as a function of the quadrupole 
deformation for othet N = 28 isotones. The binding energies of 40 Mg and 42 Si are also found 
to possess both an oblate and a prolate minimum. Just as for 44 S, we find that increasing 
the value of c pair tends to flatten their energy curves and to reduce the magnitude of the 
deformation at the minima. The calculation using c pa j r =0.55 yields, for 42 Si, a sharp oblate 
minimum at /3=-0.32, about 1 MeV below a shallow prolate minimum at (3=0.25. For 40 Mg, 
the calculation using c pa j r =0.55 yields a fairly shallow prolate minimum at /5=0.38, about 
0.5 MeV below a shallow oblate minimum at (3=-0.25. As far as we know, the deformations 
of these nuclei have not been measured. As mentioned before, a small neutron pairing gap 
persists in our calculation of the spherical nucleus 48 Ca. This gap disappears in the heavier 
iV=28 isotones, 50 Ti, 52 Cr, and 54 Fe, which are also spherical, but which possess instead 
an average proton gap. In Table I, we provide a summary of the ground state energies, 
deformations, and average gap parameters of the nuclei in the N = 28 chain, obtained 
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with Cp a i r =0.55. We note that, except for the lightest nuclei in chain, 38 Ne and 40 Mg, the 
calculations are in good agreement with the Audi-Wapstra and Moller-Nix systematics. 

In agreement with the other calculations discussed above, we find the even-even nuclei 
with A < 48 in the N = 28 isotonic chain to display the characteristics of shape coexistence. 
The nuclei vary between oblate and prolate ground state deformations and possess both 
prolate and oblate minima that are very similar in energy. The values we obtain for the 
ground state deformations are very similar to those of the other calculations. The values 
obtained for the pairing gaps are different however. The RMF+BCS calculations of Refs. 
and |5T| fixed the gap at either 75 keV or at 500 keV, values which are relatively small 



for such light nuclei. The RHB calculations of Ref. ||42|| , which are the most similar to 
our calculations in terms of method, obtained average values of the pairing gap that are 
about 40% larger than ours, but values for the deformations that are almost identical to 
ours. If we were to increase the magnitude of the pairing field in our calculations, the 
resulting deformations would be smaller than those of the RHB (and other) calculations. 
This discrepancy is an indication of the relativistic effects of the Dirac pairing field that 
were discussed in the Introduction. There, we suggested that the Dirac structure of the 
pairing field should result in its being more localized on the nuclear surface and, thus, more 
effective in limiting the deformation. As the Dirac structure of the pairing field is the 
principal difference between our calculations and the RHB ones, our results lend support to 
this argument. 



C. Deformed nuclei in the Z « 40 region 

Nuclei in the region of the Z « 40 subshell closure present interesting variations in 
deformation as they deviate from the neutron magic number N = 50 on either side of the 
stability line. Studies of their level schemes and lifetimes have shown that the ground-state 
properties of these nuclei are very sensitive to small changes in the proton and neutron 



number f|63l|64 



Strontium isotopes present two regions (N 38 and iV « 60) of very strong deforma- 
tions, as large as /3 = 0.4. This behavior can be seen in the E2+ measurements as well as in 
the isotope shifts. The rapid transition from a spherical shape near the iV=50 shell closure 
to a strongly deformed one can be explained by the reinforcement of the shell gaps in the 
single-particle levels for protons and neutrons [|]]. 

For krypton isotopes, on the other hand, the expected change in shape at iV = 60, which 
occurs for most elements in the Z = 40 region, does not appear. Instead, the isotope shifts 
show a pronounced slope change at N = 50. The fact that a strongly deformed ground state 
does not appear at N = 60 has been attributed to a reduced proton-neutron interaction at 
Z = 36. 

Theoretical studies of these nuclei have been carried out using both relativistic and non- 



relativistic frameworks. Bonche et al. |54] , |55|| studied these isotopes extensively including 



triaxial deformations in the non-relativistic Hartree-Fock with a Skyrme SIII force, which 
resulted in a linear trend for the isotope shifts. They also used the method of generator 
coordinates for Sr isotopes and predicted a transition from a spherical to a deformed shape. 
However, the isotope shifts could not be reproduced in any of those calculations. 
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Many studies using the relativistic mean field approach have been performed in this re- 
gion. Early calculations employed the parameter set NL1, which has a very large asymmetry 
energy and does not describe the isotope shifts well [20 , 21] . In Ref. pZ |, however, Lalazissis 



and Sharma successfully described the ground state properties of the Sr and Kr isotopes, by 
applying the RMF approach with the NLSH parameter set for the mean field Lagrangian and 
the BCS formalism for the pairing correlations. The deviations of the binding energies from 
the experimental data were the on order of 0.5% and the deformation parameters j3 were in 
good agreement with the (3 values extracted from B(E2) measurements. It is therefore of in- 
terest to perform calculations of the Sr and Kr isotopes using the Dirac-Hartree-Bogoliubov 
approach in those regions and compare the results with the experimental data. 

We begin by studying the dependence of the binding energy on the quadrupole defor- 
mation parameter, (3, using the method described in the previous section. In Fig. |T0, we 
present the total binding energy as a function of the deformation (3 of the nucleus 100 Sr, for 
three values of the parameter c pa i r . The most singular difference between the three curves is 
the change of the absolute minimum, from a prolate shape in the case of c pa i r = 0.51 to an 
oblate one for c pa i r = 0.54, remaining oblate for c pa i r = 0.58. The energy difference between 
the two minima is quite small and the oblate minimum becomes deeper relative to the pro- 
late one as c pa i r increases. For c pa i r = 0.51, the prolate minimum is 379 keV deeper than the 
oblate one. For c pair = 0.54, the oblate minimum is 17 keV deeper, while for c pair = 0.58, it 
is 823 keV deeper. If we take the oblate solution as the absolute minimum, the value of the 
deformation parameter agrees well with the experimental value of (3 = 0.372 ± 0.008. 

We performed calculations for possible minima of both prolate and oblate shapes of most 
of the known even nuclei in the Kr and Sr isotopic chains. We found many nuclei for which 
the two minima have very similar binding energies, an indication of shape coexistence or of 
a triaxial ground state deformation. 

The two-neutron separation energies of the even Kr and Sr isotopes, as functions of the 
mass number, are shown in Fig. O, for three values of c pa i r - To prepare this figure, we 
have selected the solution, among the multiple minima of each isotope, that has the deepest 
(absolute) minimum. We compare our results with the Moller-Nix systematics and the Audi- 
Wapstra compilation. We find the calculations to agree quite well with both of these. As in 
the previous cases, the calculated two-neutron separation energies are fairly insensitive to 
the parameter c pa i r . 

In Fig. [l^, we display the mean value of the neutron gap parameter < A n > of the even 
mass Kr and Sr isotopes, as a function of the mass number, for the same three values of 
Cpair- We find the best agreement with the Moller-Nix and Audi-Wapstra values to be given 
by the calculations using c pa i r = 0.58. The shell closure at iV = 50 is clearly visible as a zero 
in the calculated pairing gap of both isotopoic chains. 

In Table II, we present the results for the total binding energy and the deformation 
parameter f3 obtained in the present calculations, for c pa j r =0.58, together with the experi- 
mental data. The value c pa j r =0.58 furnishes the best description of the total binding energies 
and of the mean value of the neutron pairing gap. The cases in which the difference be- 
tween two minima is less than 200 keV are marked with a star and the results for the two 
deformations are presented. In general, the calculated binding energies agree quite well with 
the experimental values, although the proton-rich isotopes are slightly underbound in the 
calculations. 



27 



Our calculations, with one exception, predict that the non-spherical even nuclei in the 
two isotopic chains will have an oblate ground state deformation. This contrasts sharply with 
the Moller-Nix systematics which, with two exceptions, predicts prolate deformations for the 
ground states of the non-spherical nuclei of the two chains. Such disagreement is not too 
surprising, given the sensitivity to the interaction parameters of the binding energy versus 
deformation curve of many of these nuclei. An example of this sensitivity was displayed in 
Fig. [TT], in which the ground state of 100 Sr was seen to be prolate for c pa i r = 0.51 but oblate 
for c pa i r — 0.54 and 0.58. We note that, of the the five nuclei for which we found a second 
minimum within 200 keV of the ground state one, the deformation of the second minimum 
of four of these nuclei is in agreement with the Moller-Nix prediction. Experimentally, large 
ground state deformations of (3 ~ 0.44 have been found in 76 Sr, 74 Kr, and 76 Kr ||]. These 
values are in agreement with the Moller-Nix ones and with our result for 76 Sr. They are not 
consistent with our results for the Kr isotopes. 

The calculated isotope shifts of the even-even Kr and Sr isotopes, with respect to the 
reference nuclei 86 Kr and 88 Sr, respectively, are compared to the experimental data in Fig. [13. 
We display only the results for the deepest energy minimum of each isotope in the figure. 
In the cases in which a second minimum is close in energy to the deepest one, their isotopic 
shifts are also very similar. We find the agreement with the data for the neutron rich side 
to be very good. In addition, the positions of the kink due to the shell effect for the Kr and 
Sr isotopes are also well reproduced. However, the calculations underestimate the isotope 
shift on the proton rich side of the stability line. 



VII. SUMMARY 

We have used an extension of the Gorkov formalism for the description of pairing to 
develop a Dirac-Hartree-Fock-Bogoliubov approximation to the ground state wave function 
and energy of finite nuclei. We have applied it to spin-zero proton-proton and neutron- 
neutron pairing within the Dirac-Hartree-Bogoliubov approximation (we neglected the Fock 
term). We have retained the full Dirac structure of the pairing field that is permitted by 
the symmetries of the problem. We find the Dirac pairing structure to be dominated by a 
scalar term and the zero component of a vector term, as has also been found to be the case 
for 1 So pairing in symmetric nuclear matter. 

In our calculations, we use a zero-range approximation to the pairing interaction, which 
results in a local pairing field. We justified the zero-range approximation by arguing that 
the effective length for spatial variations of the wavefunctions, in the calculations performed 
at present, is larger than the range of the nonlocality of the pairing interaction, rendering 
the effects of the nonlocality close to negligible. 

We studied the effects of the Dirac pairing field on the properties of the even-even nuclei 
of the isotopic chains of Ca, Ni and Sn (spherical) and Kr and Sr (deformed), as well as 
the N=28 isotonic chain. We first studied the isotopic chains of the spherical nuclei in 
order to determine the sensitivity of various nuclear observables to the interaction in the 
pairing channel. We found the two-neutron separation energies, root-mean-square radii and 
isotopic shifts of the spherical nuclei to be fairly independent of the strength of the pairing 
interaction. The binding energy, on the contrary, is quite sensitive to the pairing interaction 
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strength (at least when the pairing is not identically zero due to a shell closure). Quite 
obviously, the average pairing gap also displays a strong dependence on the strength of the 
pairing interaction. We observed that an adequate choice of this strength simultaneously 
yields reasonable values for the binding energy, near the stability line, and rough agreement 
between the average pairing gap and the experimental gap, defined in terms of even-odd mass 
differences. Far from the stability line, our calculated binding energies showed increasing 
discrepancies with the experimental data, suggesting that the NL3 parameter set might 
require further adustment, at least when used in a DHB formalism. All in all, however, our 
results show good agreement with the data and with previous calculations. 

To study the deformed Kr and Sr nuclei and the deformed nuclei of the N=28 isotonic 
chain, we performed a sequence of constrained Dirac-Hartree-Bogoliubov calculations, as a 
function of the quadrupole deformation (3, for each of the nuclei in question. Each sequence 
of calculations provided us with a curve of binding energy versus deformation that permitted 
the localization of the equilibrium configurations of the nucleus. Most of the nuclei in these 
three chains present both a prolate and an oblate minimum (and, at times, a slight spherical 
minimum as well). In agreement with other calculations, we found nuclei in each of the 
chains for which the energies of the oblate and prolate minima were very similar, raising the 
possibility of shape coexistence or triaxial ground state deformations. For the most part, 
our results present reasonable agreement with the data and with previous calculations. We 
did, however, observe a discrepancy between our average pairing gaps and those of the RHB 
calculations of Ref. [[42], which we took as an indication of the relativistic suppression of 



the pairing field due to its Dirac structure. Our DHB calculations furnished quadrupole 
deformations similar to the RHB ones for average values of the pairing gap that are about 
30% smaller than theirs. 

Based on the nuclei we have studied here, we conclude that the DHB approximation can 
provide a description of the binding energies, rms radii and ground state deformations of 
finite nuclei that is at least as good as that provided by the nonrelativistic HFB approxima- 
tion or by the RHB approximation. Our claim that the DHB approximation can provide a 
more reliable description than the others, however, has not been demonstrated here. This 
will require further study of the parameters that enter the pairing interaction, in order to 
better fix their values. 

In the future, we plan to extend our calculations to a larger set of isotopic chains and to 
analyze in detail their proton- and neutron-rich tails, in order to obtain a set of parameters 
that describes these better. We plan to extend our calculations to odd nuclei by including 
a blocking term. As the blocking term varies, along with the other states, during the 
search for the Dirac-Hartree-Bogoliubov stationary point, this does not seem to be a trivial 
task. In order to facilitate it, we plan to first develop an extended BCS approximation, 
based on the self-consistent Hartree eigenstates of the mean field in conjunction with the 
relativistic pairing field. Such an approximation could also serve as a faster first stage to 
the Dirac-Hartree-Bogoliubov calculations, which are quite time-consuming at the present. 
Eventually, we also intend to include the nonlocality of the pairing interaction, in order to 
use more realistic interactions containing the relatively long-range effects of pion exchange. 
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APPENDIX A: TIME-REVERSED PAIRS OF STATES 



We begin by writing the complex conjugate of Eq. fl25p as 



V BA*(x,$)Bi ((s-n)6(x-&)+Bh* T (x,y)Bi) 7o J \BV*(y) J 

where B = 75C is the Dirac matrix part of the time-reversal operator. Analyzing the 
Hamiltonian term, we see that 

B h*(x, y)B ] = B {ia* ■ V + l3*M)B^5(x — y)+B /?*£*(£ , y)B ] (A2) 
= (-ia- V + f3M)5(x -y) + B P*X*(x,y)B* . 

The self-energy term, in turn, can be put into the form 

B^(x,y)B^ =8(x-ff)J2 B ^U S ) Bi [ d 3 zDf(x-z) £ [z)B^ B^T^B^ BU;{z)B^ 

3 £ 7 <0 

-J2 B ^* a (x)B^D^(x-y) BU;(x)U^(ff)B^B 7 * T^(y)B^. (A3) 

3 e 7 <0 

If we now look closely at the vertex functions we are considering here, we find 

for the a meson: B%\B^ . . . B^IB^ — 70 1 . . . 7 1, 

for the u meson: B%%B^g^ . . . B^*j*B^ = 7o7 M 9^ ■ ■ • 7o7^, 

for the p meson: B~/*j*t* : B ] '5 ij . . . B%1* v t*B^ = 707^ g^8 ij . . . 7o7^j, 

for the photon: B^fi-^B^g^ . . . B^fi-%^ = i±±^ 7o7/ ^ . . . ^ 7o7 ,. In 
short, the products of the pairs of vertices remain unchanged by the transformation. We 
thus have for the transformed self energy, 

BP* E*(af,0flt = 5(x - y) J2loT ja (x) f d 3 zDf(x-z) I £ (z)B^ T 3 p(z)BU;(z) I (A4) 

3 J \e-,<0 J 

-Y,7oT ja (S)D^(x-g) ( J2 BU*(x)U^(y)B^ j T jl3 (y) , 

3 \e T <0 J 

in which the only differences from the original expression for the self energy are the prop- 
agator contributions, written in terms of the solutions of the Dirac-Gorkov equation, that 
enter the expression. 

We can perform the same analysis on the pairing term, for which we find a similar result, 
BA'&jftB* = -£r ia (f)I}f ( £ BU;(x)V?(y)B^A AT%(y)A^ , (A5) 

j \e 7 <0 / 

in which the only difference between the original expression and the transformed one is again 
the propagator's contribution to the expression. Analyses of the other terms, B h^(x, y)B^ 
and B A*(x, y)B\ yield similar results. 
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We can thus show that the transformed wave vector ( BU ^ | satisfies a similar equa- 

\BV*(x)J 

tion 

d3 ( lo((e + ^)5(x-y)-h(x,y)) A(x,y) \fBU*(y)\ Q 

V \ A(x, y) ((e - »)6(x - y) + h T (x, y)) 7 o J \ BV*(y) ) ' 

to that satisfied by the wavevector • Each is a solution of an equation with the 

y V W J 

same energy eigenvalue but with orbital quantum numbers and mean fields that have the 
time-reversed Dirac structure of the other. If each member of these pairs of states possessing 
time-reversed Dirac structure are equally occupied, the sums over states that enter into the 
definition of the self-energy and pairing fields will be invariant under the transformation 
that temporally inverts the Dirac structure. That is, we then have that 

E U,(x)U^ = J2 BU;(x)UT(y)B* (A7) 

£ 7 <0 £ 7 <0 

and 

E U^V^y) = E BU;{x)V^{y)Bl (A8) 

£ 7 <0 £ 7 <0 

In this case, both of the wavevectors ( U [ X ] ] and ( BU ^f] ) satisfy the same Dirac-Gorkov 

V v ( x ) J V BV *( X > J 
equation with the same energy eigenvalue. Then either both or neither of them will enter 

the propagator's contributions to the self-consistency equations and these contributions will 

indeed be invariant under time-inversion of their Dirac structure. We assume this to be the 

case. 
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FIG. 1. Two-neutron separation energy of the even isotopes of Sn, as a function of the mass 
number A, for three values of the parameter c pa i r . The values obtained from the compilation of 
experimental masses of Audi and Wapstra |57j (solid circles) and the Moller-Nix systematics [56] 
(open circles) are also shown. 
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FIG. 2. Mean value of the neutron gap parameter of the even isotopes of Sn, as a function of 
the mass number A, for three values of the parameter c pa i r (lines). The odd-even mass differences 
obtained from the compilation of experimental masses of Audi and Wapstra |57j (solid circles) and 
from the Moller-Nix systematics [^j (open circles) are also plotted. 
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FIG. 4. Difference between the calculated binding energy and the Moller-Nix binding energy 
[56] for the even isotopes of Sn, as a function of the mass number A, for three values of the 
parameter c pa i r (lines). The difference between the calculated binding energy and the Audi-Wapstra 
experimental binding energy [57], for c pa j r =0.500, is also plotted (circles). 
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FIG. 5. Deviation of the root mean square neutron radius from the systematic value roiV 1 / 3 , 
for the even isotopes of Sn, as a function of the mass number A, for three values of the parameter 
c P air- The reduced radius, ro, was adjusted to the root mean square neutron radius of 132 Sn. 



pair 




FIG. 6. Mean value of the neutron gap parameter of the even isotopes of Ni, as a function of 
the mass number A, for three values of the parameter c pa i r (lines). The odd-even mass differences 
obtained from the compilation of experimental masses of Audi and Wapstra |57j (solid circles) and 
from the Moller-Nix systematics [^!J (open circles) are also plotted. 
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FIG. 7. Average value of the gap parameter of the even isotopes of S, as a function of the mass 
number A, for two values of the parameter c pa i r (lines). The odd-even mass differences obtained 
from the compilation of experimental masses of Audi and Wapstra |57[ (solid circles) and from the 



Moller-Nix systematics [56| (open circles) are also plotted. 
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FIG. 8. Binding energy of 44 S, as a function of the quadrupole deformation f3, for two values 
of the parameter c pa i r . 



40 





41 



30r 





% 80 90 100 % 80 90 100 

A ( A 

FIG. 11. Two neutron separation energy of the even isotopes of Kr and Sr, as a function of the 

mass number A, for three values of the parameter c pa i r - The values obtained from the compilation 

of experimental masses of Audi and Wapstra |57j (solid circles) and the Moller-Nix systematics 

f)6| (open circles) are also shown. 
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FIG. 12. Mean value of the neutron gap parameter of the even isotopes of Kr and Sr, as a 
function of the mass number A, for three values of the parameter Op a i r (lines). The odd-even mass 
differences obtained from the compilation of experimental masses of Audi and Wapstra (solid 
circles) and from the Moller-Nix systematics p6| (open circles) are also plotted. 
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FIG. 13. Isotopic shifts of the even isotopes of Kr and Sr, as a function of the mass number 

A, for c pa j r =0.58 (lines). The experimental values are displayed as points with error bars. 
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